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L. SYNCHRONOUS MACHINE MODEL DEVELOPMENT 


The primary means of shipboard electric power generation is with synchronous 
generators driven by steam turbines. The electrical and electromechanical behavior of the 
synchronous generator can be accurately predicted from the equations for a three-phase 
salient-pole synchronous machine. 

The rotor of the synchronous generator is equipped with a field winding and one or 
more damper windings. The rotor windings usually all have different electrical 
characteristics. Besides the differing electrical characteristics, the salient-pole machine 
exhibits a magnetically asymmetrical geometry due to the varying air-gap path for the 
magnetic flux. This magnetic asymmetry creates time-varying self- and mutual-inductance 
terms in the voltage equations for the stator variables. By transforming the equations into 
a reference frame which rotates at the same angular velocity as the rotor, the magnetic flux 
path remains constant, thereby eliminating the time-varying inductance terms. 

This chapter develops the voltage and electromagnetic torque equations for the 
synchronous generator in terms of machine variables and then converts these equations to 
the rotor reference frame via Park's transformation [1]. The developed model assumes a 
magnetically linear system thereby neglecting the effects of saturation. Saturation could 
be accounted for with the addition of another simulation segment for the model, however 


saturation data for the machines studied was unavailable. 


A. VOLTAGE EQUATIONS IN MACHINE VARIABLES 

Voltage equations for the synchronous generator will be developed from a simple 
two-pole, three-phase, wye-connected, salient-pole synchronous machine shown in Figure 
1. The number of poles of the particular machine can be easily changed as this 1s one of 
the variables of the resulting torque equations. The stator windings are assumed identical 
with sinusoidally distributed windings, displaced 120°, with N, equivalent turns and 
equivalent resistance r, (subscript s denoting stator variables). A field winding and three 
damper windings (fd, kd, kg1, and kq2 respectively) are installed on the rotor. All 
windings are assumed sinusoidal with applicable N equivalent turns and resistance r 
(subscript denoting particular circuit). The quadrature axis (g axis) and direct axis (d axis) 
are used with reference to the magnetic axis of the rotor circuits. Voltages of the rotor 
damper windings (v,,;, ¥,,,, and v,,,) are for monitoring purposes only, the windings 
themselves are physically shorted. 


Voltage equations in machine variables [2] may be expressed in matrix form as 


V abes = — Vd ahes + PMabes (1) 


V oar =Vriga + Doar (2) 


anes _ L, Le hela 
porte ole ° 


where 





Figure 1. 2-Pole, 3-Phase Synchronous Machine. 


and 


L, +L, —L,cos20, - sh, — £, cos 2(0, a) - sh, — L, cos2(0, ee 
L, = -sh, — L, cos2(0, - a L, + L, — L,cos2(@, — =) = sh, ~L,cos2(0, +1) 
- xh, — L, cos2(0, ad = sh, -~1,cosX6,+n) L,+L,—L,cos2(0, + =) 
(4) 


L,,, C089, L492 C0S8, Ly,sin 8, Ly4 sin 9, 


L,= L008 (8, -=) Lg 0080, -—) Lj, sin @,-=) Lagsin @,-—) (5) 


Digg COS 0,422) Lys 008 (0, +=) Ly Sin (0,422) Los sin(0, +72) 
3 3 3 3m 
| Ea oa Hae Lig tha 2 0 0 
Ta Lg 1hq2 Lg? 1 Lng 2 0 ° (6) 
’ 0 0 Liye BE Lf Lexa 
0 0 L paz Loy ae 


In (1) and (2), the time differential operator is p. In (4) and (6), the subscript / denotes a 
leakage inductance term. The subscripts skgl, skq2, sfd, and skd denote mutual 
inductances between stator and rotor windings. 


The magnetizing inductances are defined as 


Lng = = (Ly ~ Ly) (7) 


li = S(L, + Ly) (8) 


where the term incorporating L, refers to the inductance term which is independent of 


rotor position and the term with Lp is the inductance term which is position dependent. 


The rotor variables are then referred to the stator windings via the appropriate turns 


ratio giving the following equations 
bow oe a (9) 
, a \T t op; 
os 3 (Le ) L. J oar 


Lg COSO, Lng COS®, L,sin8, L,sin@. 
Li =| 4,, cox8, -=) L.._ coX8, -) L, sin(8. = I, sin(8, -) (10) 


|b ae oe ie 0 0 


L’ = Des. orp ae . 0 0 (11) 
0 O Latha Ln 
0 0 ee Lig + Lng 


where the primed variables (i.e. Li. ) refer the appropriate rotor variable to the stator 
windings. 
The voltage equations expressed in terms of machine variables referred to the stator 


windings are 
V r, as pL, pL, —] 
a = Z i T , L’ a (12) 
V adr 3 P( 7) I, + Pp r Nod 


Note that the voltage equations (12) are valid for positive direction of stator current 


assumed out of the stator terminals. If motor action is assumed, the sign preceding i,,.. 


must be changed. 


B. TORQUE EQUATION IN MACHINE VARIABLES 
The field energy stored in the coupling field of a synchronous machine can be 


expressed as 


ae : a 
W, = 5 ates YP (L, — Li Dijge — Gate)? Li 


sr’ qdr 
i (13) 
+5 Gite (Li - LAD, 
where 
Li, =diag| Lig, Ligr Lye Lia] (14) 


The electromagnetic torque is found by taking the partial derivative of the coenergy with 
respect to the rotor angular displacement 8,. The magnetic system is assumed linear, 
therefore the coenergy equals the field energy, and upon taking the partial derivative, the 


electrical torque is 


mad 


a0 [L, a Lr. Vigses + (iases) 


te Ph ae 2 uKe| (15) 


where P is the number of poles of the machine. The torque is positive for generator action 
with the positive direction of stator current assumed out of the stator terminals. 


The torque and rotor speed are related by 


T= -() po, +T, (16) 


where J is the inertia of the machine in kilogram meters? (kg-m2) or Joule seconds? (J:s2). 


The input torque T; is positive for a torque applied to the shaft of the synchronous 


machine. 


C. VOLTAGE EQUATIONS IN THE ROTOR REFERENCE-FRAME 

The machine voltage equations (1) and (2) are seen to contain flux-linkage terms 
which are dependent on rotor position 8,. The analysis of the machine behavior is greatly 
simplified by transforming these time-varying inductances to constant inductance terms. 
This can be realized if the stator variables are transformed to the rotor reference frame. 
R.H. Park was the first to incorporate a change of variables in the analysis of synchronous 
machines [1]. 

A change of variables which formulates a transformation of the 3-phase variables of 
the stationary circuit elements (thus the need to refer the rotor variables to the stator) to 


the rotor reference frame may be expressed as 


0s < K faves (17) 
where 
te 
£405 7 Sas (18) 
hos 


Eaves a Ss (19) 


and 


cos6_ cos(8, -) cos(8, +=) 


i Se ge 
K' = ; sinQ, sin(0_ 3 ) sn(@.+ : ) (20) 
m u 
2 2 2 
6, = ['o,€) &+0, 0 (21) 


The term K* denotes the transformation matrix of stator variables (the subscript-s) to the 


rotor reference frame (the superscript-r). It can be shown that the voltage equations (1) 


and (2) are transformed into 


V d0s = =z 0: av. Nias + PNia0s (22) 
vi, = i, + Pace (23) 
where 
Na 
QW) = |e (24) 
0 


The flux-linkages can be expressed in the rotor reference frame by transforming (9) to 


obtain 


r ry) ry 7? 
N ad0s a eee KL, mG: 
Mie] |F@ KL 


r 


Expanding the voltage equations given by (22) and (23) gives 


vi, = Td), +O, + PN, 
Vig = Td, —O,N,, + PN, 
Vos = Tying + Phos 

Voor = Negilign + Phi 

Vega = Negatiga + PMig 

Vig =Nalgg + PM, 


Via =Teabeg + PMea 
and the expanded flux linkage terms become 


Nog =— yin + Lang (is, His + iff) 
Nas = — Lights + Lng (tas + tea FH) 

Mos = Lyles 

r’ i= Lugrtegs Bs less Gre + ipa 7 lpg?) 
Na Legale + Lag (“Ee +t iff) 
ge Lith + Ls (is, tig + 


NM a= Lnaieg + Lina (tg + fg + tg 


(25) 


(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(32) 


(33) 
(34) 
(35) 
(36) 
(37) 
(38) 
(39) 


Expressing the voltage and flux linkage equations in terms of reactances rather than 
inductances gives the voltage equations in matrix format (40), where the expression for 


the applied field voltage is introduced 





Tr oT Xx 
eft, =v nd (41) 


f 


fd 


and the leakage and magnetic reactances are combined as 


2 pnt Kp (42) 
X,=X,+X_, (43) 
x, = Xe (44) 
Nie = ae (45) 
X)) = eee (46) 
X pf Se (47) 


The zero sequence terms have been omitted since balanced conditions are assumed. 

The model being developed will be connected to a non-infinite bus. This implies that 
the voltage and frequency do indeed vary with load, based on the droop characteristics of 
the particular regulator. With this in mind, it is necessary to solve for currents based on the 
bus voltages, therefore the currents are selected as independent variables. Rewriting (40) 


with the current and current derivative terms separated. 
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(48) 


The state equation format of (48) is used for the simulation, so the current derivatives 


must be solved for. As an example 


v= [A |i +[Blpi 


pi=([B] ‘{v —[A]i} 


and, with this format, the current derivatives are solved as shown in (50) 


12 


(49) 


(0S) 





















































0 
0 
0 
vas 
q 
bu © 
X77 - 
0 
0 0 
0 0 
Ig 1M 
ory sa) 
1@ 1M 
shes [oy 
0 0 
io a @ 
“ide “aN. 











IM 





oO oO & & 











13 


Denoting the inverse matrix elements as [B] * = [b] = 
O61 


the desired differential current equations are obtained 


7 = OylU,, FEL Aa jr Mod jor o Sona) 


Ls eet OR ee r 
W, b b 


, Ate 
a5 Oi3l Veer = Tegal 
ie 
+ B,4[Uy2 — Nigalige | 
pift, = byLuf, + ri, + 0, Bi, - @, if - @, Am i) 
@, b @, 
+ OsslUpa — Tetteg | 


Oey 5 Tegolice | 


ald, + Ay ar And irr gy And. jr 


W, ®, ' Os 


+ OlUy — Tigtea | 


+ OglUy, — Nala | 


7 ae r Tr 
Dk, si byo[U3, a Teas -@ 


r qs 


ar by <Levfa — Xda] 
+ Bog lUia — Tedea | 
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(51) 


(52) 


(53) 


(54) 


Xing et Sng ist 
+ @ a Oe bigot 
w w 


XxX, . 
Dlg = Osolvg, + Tg, — O, — 1 
WO, b b 


ap ble ry = nda! (55) 


+ bsg [Up — Tealea | 


2 cal r St ° “or 
Dig = OgelUz, + Tlgs — Saat Ligs + O 


r 


+ Ogle via — Xndfa | (56) 
Son eae | 
Note the similar form of the g- and the d-axis integral current equations. These will 
greatly simplify the SIMULINK [3] model. 


D. TORQUE EQUATION IN SUBSTITUTE VARIABLES 


The torque equation (15), transformed into rotor reference frame variables is 


eles i: ; Oe eee te tee 5 5 
f= (SMS na tas TY Uo) Be teat Bod (Se ees pe eye io?) 


and writing in terms of reactance 








(58) 
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gives 


3 | x “7 of ofr bat X bd °/7T °/r “7 
Oho (=i, iy oe tet) Lal (59) 


1, Sul 


The rotor speed is obtained from (16) as 
w, = (-——-T-T) (60) 
Dp r VJ e I 


E. ROTOR ANGLE 
The rotor angle of a synchronous machine is a useful parameter to monitor the torque 
production of the machine. It relates the angular displacement of the rotor and stator 


magneto motive force axes. In equation form 


5=6, -6, 
| (61) 
= flo, (8) ~ w,,(§)]d& + 6,(0) — 6,, (0) 
where 6 is the rotor angle (in radians), and € is a dummy variable of integration. 


Equation (61) can be implemented by solving equation (61) in terms of the available 
parameters. The rotor position, 6, is available by integrating the rotor speed. However, 


the voltage angle 6,, must be solved for as follows 


0, = [o,(&) & +6,,(0) (62) 
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and 


< 
lI 


as = Vg, COSO, + v,, sin 0, 


ey 
ne + v,, cos[@, + tan Gop)! 


qs 


= a Ya + Vi cost | @, dt + 8,(0) + tan”! (—“)] 
VY 


qs 


Se 


) 


a 


(63) 
therefore, from (62) 
@., = [o,(&) dé +0,,(0) 
= | @, dt + 0,(0) + tan"! (4) (64) 
@, # 
giving 
=@ -6, = tan}(-“) (65) 
y 


qs 


where a positive rotor angle is assumed for generator action. 
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i. SPEED AND VOLTAGE REGULATORS 


For the machine model to accurately simulate a shipboard power system it is 
necessary to include controllers for both the field voltage and the prime mover torque. 
The two principal controllers of a turbine generator set are the automatic voltage regulator 
and the speed governor. The automatic voltage regulator model consists of voltage 
sensing equipment, comparators and amplifiers controlling a synchronous machine which 
can be generating or motoring. The speed governor may be considered to have simular 


equipment but in addition it is necessary to take the turbine into account. 


A. SPEED REGULATOR 
1. Model Development 
The block diagram of a steam turbine and associated governor for a typical 
naval ship's service turbine generator (SSTG) is shown in Figure 2 [4]. 
The reference speed input is the no-load frequency selected by the operator of 


the electric plant. The gain K is used to adjust the speed droop of the regulator. T, is the 


time constant associated with the servo and steam valve system andT, is the time constant 


associated with the conversion of steam to torque (T7,,, ). The following values were used 


for this study 
15, 0.1 sec 
oe 0.06 sec 
K 150 


From this block diagram (Figure 2), a SIMULINK model was developed and is shown in 
Figure 3. 
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Figure 2. Block Diagram, Steam Turbine & Governor. 
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Figure 3. SIMULINK Speed Regulator. 


As can be seen, the speed regulator is designed to compare the reference speed with the 
speed being produced by the synchronous machine. This speed difference determines the 
prime mover torque being applied to the generator. When the speed difference is no 
longer changing, the prime mover torque is equal to the electrical torque developed by the 
generator (load dependent). The final machine speed (and therefore the machine output 
frequency) will be at some value below the no-load set-point frequency. This can be 
demonstrated by referring to the ‘house’ diagram shown in Figure 4. 

As the output power of the machine is increased, the operating frequency of the 
machine decreases (based on the slope or droop of the particular regulator) the operating 
frequency can be restored to a desired frequency (i.e., 60 Hz) by the operator adjusting the 
no-load setpoint of the regulator. 

2. Regulator Stability 

As was mentioned in the previous section, the speed droop characteristics of the 
regulator are determined by the gain K of the regulator. The percent speed regulation is 
defined by 

-O, 


0) 
% Speed Regulation = —~—* 100% (62) 
w 
fl 


As the gain K is raised, the percent speed regulation approaches 100%, at which 
point the regulator and therefore the regulated plant will become unstable. This stability 
point was determined by examining the root locus of the model of Figure 3. The 
synchronous machine model was obtained by linearizing the machines current block 
around its stable operating point for rated conditions by using the MATLAB [5] function 
linmod. The linearized plant model returned by linmod is in linearized state space format 


which is then changed to transfer function format with the function ss2¢f. The transfer 


21 


IoMOg 
suneisdg 


(MY) IOMOd 0 


Aouonboly 
sure1sdO 


Aouonboljy 


Figure 4. House Diagram. 
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function format of the plant, speed block, and speed regulator is then reduced to a single 
open loop to study the roots of the system as the regulator gain K is varied. This 
complete process is demonstrated in Figures 5, 5a, 5b, and 5c. The resulting root locus 
plot is shown in Figure 6. The operating constraint was determined for the first root 


meeting the imaginary axis for a given K. The limit of K was determined to be 153. 


23 


Figure 5. Synchronous Machine Model. 
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Figure 5a. Current Block. 
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Figure 5b. Linearized Current Block, Speed Regulator, and Speed Block. 
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Figure 5c. Transfer Function Development, Speed Regulator. 
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Figure 6. Root Locus Plot--Synchronous Machine and Speed Regulator. 
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B. VOLTAGE REGULATOR 
The voltage regulator or excitation system used in this study is shown in block 


diagram form in Figure 7 [6]. The parameters are defined as follows 


t, and K, Regulator amplifier gain and time constant 
t, and K,; Saturating exciter gain and time constant 
cr) Ceeand i. Feedback loop gain and time constants 


The SIMULINK model developed from the block diagram is shown in Figure 8. The 
regulator compares the bus voltage with the specified reference voltage. As can be seen in 
the Simulink model, this comparison is done in the gd0 reference frame. 


The reference voltage is computed as follows 
Vieference a y Vy SE Vg ; (63) 
This equates to the desired line-to-neutral voltage (peak) as shown below 


os 
V ad0s a ARE (64) 


cos@, cos(0, -) c0s(0, +) 


Vy Va 
vi = sin8, sin(8, ails sin(6. pe v, (65) 
3 3 3 
7 a 1” Ly. 
2 2 2 = ui 
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Figure 7. Block Diagram, Voltage Regulator. 
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Figure 8. SIMULINK Voltage Regulator. 
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where 


v, =v, cos(8,,) 


v, =v,, cos(8,, >) 


v. =v, cos(8,, + 


and 


8, = ['0, © &d+0,(0) 
8, = |, @,6) Ed+0,, (0) 


and solving for Vg and vg 


V, ss cos, cos8,, +cos(0, Le), eet) Be a). +] 
3 3 3 3 3 
Vy ==, sin 8 cos8,, +sin(O. - F008 -Zy+sin(@®, +Zreos, 2) 


and using the trigonometric identities 


COS XCOS y+COSs -) Cos ( y+ COS 4) cos ( y+) == cos(x— y) 


sin x cos y+sin oS) COs ( y-)+sin (1+ D)eos ( y+) == sin (x-y) 
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(66) 


(67) 


(68) 


(69) 


(70) 


(71) 





gives 


v, ==, =cos(0, Ue) 
=v,,cos[(@, —@,)t+8, (0)—8, (0)] 
2.. 3m 
Vs Saya —§.,) 


=v, sin[((@,-@,)t+8,(0)—8, (0)) 
At synchronous speed ®, =@,, and from (61) 6=6, —8,,, therefore 


V,=V,,C0S(d) 


v, =v,,sin (6) 


Vreference = vz [cos” (5) + sin? (5)] = y,, 


and 
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(72) 


(73) 


(74) 
(75) 


(76) 


Ht. MACHINE LOAD 


To accurately study the machine's performance it is necessary to be able to load the 
given machine in a realistic manner. The inductive nature of transmission lines and 
induction motors dictate a need for a resistive-inductive (RL) loading of the machine. 

This section develops the RL-load used for the synchronous generator and then the 
SIMULINK implementation of this loading. 


A. PASSIVE RL-LOAD 

The load used for this study was comprised of a passive resistive-inductive load. 
This loading is typical of ac motors, transmission lines, electric heating, and lighting. The 
load was modeled in abc-variables in order to eliminate the need for additional 
transformations (abc — qd0). The simulation uses voltage as the input and currents as the 
outputs. This requires the phase voltages to be made up from the network equations. The 
loads also use voltage as the inputs and current as the outputs. Based on this topology it 
is advantageous to assume a Y-connected set of large resistors are connected across the 
terminals of the synchronous generator [7] and the phase voltages at the terminals of the 


generator are computed from 


v,=—Rictiz) (85) 
v,=— Rio ti) (86) 
v= — Rii.gtix) (87) 


where i,¢, i,g, and i,, are the phase currents flowing into the generator from the load and 


it> 44,, and i,, are the phase currents flowing into the load from the generator. These 
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currents are computed first as i,, i,, and i, and then transformed to i,, i,, and i,. This 
is represented in Figure 9. 


This load modeling is quite accurate as the resistance value R is made large (20-50 


pu). As the resistance is made large, the current through this branch becomes negligible. 


B. SIMULINK LOAD MODEL 
The passive RL-load was modeled in SIMULINK by starting with the voltage 
equations (only developed for phase a, others similar) from Figure 10. 


vy =e, +L, pi, +r.i, (88) 
v,=r,+L, pi, (89) 
Vv, = Ri, sh: (90) 


The differential equations can be implemented as shown in Figure 11. 

As can be seen, the terminal voltage v, is dependent on the output current i,. This 
voltage dependency is the difference between a finite and an infinite bus model of electric 
machines. An infinite bus analysis assumes no dependency between y, and i,. 

An attempt was made to make the load purely resistive. SIMULINK is unable to 
deal with this type of load however since this would incorporate direct feed-through 
between input and output of the load. This forms an algebraic loop and cannot be 
analyzed. The inductance of the transmission line overcomes this limitation. The purely 
resisitve load, connected by a transmission line (with inherent inductance) to the generator 
bus adds sufficient inductance to warrant modeling. This inductance then incorporates a 


time-delay between input and output due to the integral form of the load equations. 
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Figure 9. Generator Terminal Voltage Model. 
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Figure 10. Passive RL-Load Model. 
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Figure 11. Differential Equation Implementation of Voltage Equations. 
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IV. MACHINE MODEL PERFORMANCE 


The synchronous machines performance is examined in this section. The machines 
parameters are given and the machines response to both a load change and a fault 


condition are studied. 


A. MACHINE PARAMETERS 

The machine being simulated is a 4-pole, 60 Hz, 440 V (line-to-line), 500 kW, 0.85 
power factor synchronous generator. The per unit parameters on the machine base (Z base 
= ().1291 ©2) are as follows 

stator resistance, r; = 0.00031 

rotor field resistance, rq = 0.00226 

q-axis damper winding resistance, rpg = 0.0224 

d-axis damper winding resistance, rp = 0.00673 

stator leakage reactance, Xj, = 0.0632 

d-axis magnetizing reactance, X;,,j = 1.4648 

q-axis magnetizing reactance, X,,. = 0.7668 

field leakage reactance, Xj¢q = 0.1545 

q-axis damper winding leakage reactance, X]¢g = 0.0771 

d-axis damper winding leakage reactance, X)¢q = 0.3676 

rotor inertia, J = 14.526 
These machine parameters are taken from the Machinery Systems Laboratory (MSL) test 
facility for a typical shipboard synchronous generator [8]. The MSL test facility is within 
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the Machinery Systems Engineering Branch at David Taylor Naval Ship Research and 


Development Center. 


B. SIMULATION: LOAD CHANGE 


The machine was initially loaded at an approximate power factor of 0.8, with the load 
calculated as follows 


Machine data: 440-V 1.5 MVA. From Figure 12 





\i| = 152e6 1 _ , ogo u4 (91) 
3. 254 

T = 1.969 ei 8) KA (92) 

= 254 e!° '36.87 ° 

Z = ——_ = 129 e?"" mQ (93) 
1. 969677" kA 

r, = 103.2 mQ X, = 77.42 mQ (94) 


With the machine loaded as above, the expected steady-state power delivered to the 
load will be approximately 1.2 MW. This does not however account for the losses in the 
stator windings. The machine then underwent a 50 % step load change by doubling the 
values of r and xX L: This was done in SIMULINK by the load-change block shown in 
Figure 13. The plots of the response are shown in Figures 14 through 21. The data 


obtained from the simulation was as follows: 


FULLY LOADED 50% LOADED 
MODEL MODEL 
Load power, MW 1.1754 0.51386 
Machine power, MW 1.1792 0.5151 
Frequency, Hz 59°99 64.32 
Power Factor 0.7992 0.7832 
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Figure 12. Initial Load Parameter Determination Circuit. 
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Figure 13. SIMULINK Load Change Block. 
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vq, 00% load decrease 
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Figure 14. 50 Percent Load Decrease, vq. 
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vd, 50% load decrease 
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Figure 15. 50 Percent Load Decrease, vd. 
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Figure 16. 50 Percent Load Decrease, va. 
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Figure 17. 50 Percent Load Decrease, field voltage. 
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Figure 18. 50 Percent Load Decrease, ta. 
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te, 50% load decrease 
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Figure 19. 50 Percent Load Decrease, fe. 
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del, 50% load decrease 
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Figure 20. 50 Percent Load Decrease, rotor angle (del). 
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wr, 50% load decrease 
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Figure 21. 50 Percent Load Decrease, rotor speed (@,). 
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C. SIMULATION: THREE-PHASE FAULT 

The stability of a synchronous machine throughout a power system following a fault 
is of major concern. A 3-phase fault or short-circuit rarely occurs and a 3-phase fault at 
the machine terminals is even more uncommon; however, it is still instructive to 
demonstrate the ability to insert a fault for a chosen time interval and study the machines 
dynamic performance. 

The machine was initially loaded at a power factor of 0.8 as in part B. A 3-phase 
fault was inserted and was left in place for 7 cycles. The fault impedance (representative 
of the bus work) is 0.001291 + 30.1291 Q. If the fault had been allowed to remain on the 
system slightly longer the machine would have become unstable in terms of not having 
returned to synchronous speed after removal of the fault or the machine would have 
slipped poles prior to returning to synchronous speed. It was determined that a seven- 
cycle fault is representative of a typical breaker reclosure timing cycle which the machine 
could have been subjected to. 

The three-phase fault was performed in SIMULINK by the fault block shown in 
Figure 22. The fault implemented by using a step function block which stepped in the 
fault impedance at a particular time and then another step function block stepped in the 
normal load at the end of the desired fault tme. The plots of the response are shown in 


Figures 23 through 30. 
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Figure 22. SIMULINK Three-Phase Fault Block. 
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vq, 7-cycle 3-phase fault 
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Figure 23. 7-Cycle, 3-Phase Fault, vq. 
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vd, 7-cycle 3-phase fault 
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Figure 24. 7-Cycle, 3-Phase Fault, vd. 
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va, /-cycle 3-phase fault 
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Figure 25. 7-Cycle, 3-Phase Fault, va. 
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field, 7-cycle 3-phase fault 
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Figure 26. 7-Cycle, 3-Phase Fault, field voltage. 
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la, 7-cycle 3-phase fault 
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Figure 27. 7-Cycle, 3-Phase Fault, ia. 
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eae te, 7-cycle 3-phase fault 
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Figure 28. 7-Cycle, 3-Phase Fault, ze. 
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del, 7-cycle 3-phase fault 
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Figure 29. 7-Cycle, 3-Phase Fault, rotor angle (del). 
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wr, /-cycle 3-phase fault 
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Figure 30. 7-Cycle, 3-Phase Fault, rotor speed (@,). 
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ve CONCLUSIONS 


A. ADVANTAGES 
The advantages to studying naval power system components in a graphical 
environment are many. To be able to see, in real time, any line or load parameter (i.e. 
current or voltage) makes the study of system component interaction very straightforward. 
The concept of "closing the loop" to make the input and output of the machine 
dependent on one another, as opposed to the infinite bus concept, is easily implemented 
with the graphical tools available. Navy power systems definitely require this non-infinite 


bus modeling due to their capacity difference as compared with commercial power grids. 


B. ANALYSIS OF SIMULINK 

Matlab's SIMULINK is an excellent tool for the type of studies conducted in this 
thesis. Model building and integration into the "power system" was easily done. 

One of the difficulties with SUMULINK included the algebraic loop problem 
necessitating an integrator in the path to provide a time delay element. In this particular 
study this was not a problem due to the resulting differential equations, however, a 
possible purely resistive load would be troublesome. 

On the issue of speed of the simulations, naturally one is always looking for a faster 
simulation. I began my research on a 486-33MHz machine, then graduated to a SPARC-2 
machine and realized a speed difference of about two times. I then moved to a SPARC-10 
and achieved a speed increase of approximately four times over the SPARC-2. Typical 
times to run a simulation of approximately 5-seconds was about 4-hours on the SPARC- 
10. In comparing this to a non-graphical simulation tool such as Advanced Continuous 


Simulation Language (ACSL) [9], a comparable model runs approximately in real time (60 
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seconds of run time data in approximately 60 seconds). The disadvantage of ACSL is, of 
course, the inability of the user to observe system output in real time and the non-modular 


modeling of systems. 


oe FUTURE WORK 

The modeling of additional load components including induction machines, another 
synchronous generator sharing the load, parallel RL-loading, and adjustable speed drives 
could all be performed in SIMULINK. As the number of integrators in the complete 
model increases, however, the run time will also increase. In addition, the existing RL- 
load can be examined under such conditions as critical clearing time, first swing transient 
stability limit, and various fault conditions such as single line to ground faults, line-to-line 
faults, etc. 

In addition, voltage regulator stability can be studied to find gain limits for this 
particular regulator. Different regulator designs could also be examined. 

Mathworks has also advertised a software package which will take existing 
SIMULINK models and compile them; this will supposedly increase simulation speed 


considerably. This product sounds very encouraging. 
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